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Abstract 

Background: The etiology of Parkinson's disease (PD) remains unclear, and environmental risk-factors such as occupation 
have attracted interest. 

Objective: The goal was to investigate occupational complexity in relation to PD. 

Methods:\Ne conducted a population-based cohort study based on the Swedish Twin Registry that included 28,778 twins 
born between 1886 and 1950. We identified 433 PD cases during the study period. Data on occupation were collected from 
either the 1970 or 1980 Swedish census, and occupational complexity was assessed via a job exposure matrix. Cox 
proportional hazard regression analyses with age as the underlying time scale were used to assess PD risk as a function of 
the three domains of occupational complexity: data, people, and things. Sex and smoking were included as covariates. 
Analyses stratified by twin pair were conducted to test for confounding by familial factors. 

Results: High occupational complexity with data and people was associated with increased risk overall (Hazard Ratio 
[HR] = 1.07, 95% confidence interval [CI] 1.02-1.14, and HR=1.10, 95% CI 1.01-1.21, respectively), and in men (HR=1.08, 
95% CI 1.01-1.16, and HR= 1.15, 95% CI 1.03-1.28, respectively). Complexity with things was not associated with risk of PD. 
When the analyses were stratified by twin pair, the HRs for occupational complexity with data and people were attenuated 
in men. 



Conclusions: High complexity of work with data and people is related to increased risk of PD, particularly in men. The 
attenuation of risk observed in the twin pair-stratified analyses suggests that the association may partly be explained by 
familial factors, such as inherited traits contributing to occupational selection or other factors shared by twins. 
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Introduction 

Parkinson's disease (PD) is the second most common neurode- 
generative disorder in adults aged 60 and older [1]. With the 
continued growth of the older adult population, PD prevalence in 
people over age 50 is expected to double from between 4. 1 and 4.6 
million in 2005 to between 8.7 and 9.3 million by 2030 in 15 of the 
world's most populous nations [2]. Designing effective prevention 
and treatment options is directly linked to understanding risk 
factors for PD, making risk factor research essential. Etiology of 
PD is likely influenced by both environmental and genetic factors. 
Data from several studies suggest that environmental factors may 
be more important than genetic factors in risk of PD [3]. 
Moreover, environmental exposures are often potentially modifi- 
able and can be especially useful in disease prevention. Thus, it is 



unsurprising that numerous aspects of the environment have 
already been investigated with evidence accumulating that 
occupation-related exposures may be important [4—6]. 

Notably, while Alzheimer's disease research unequivocally 
indicates occupations reflecting low socioeconomic status [7] or 
low complexity of work [8] as risk factors, the picture is less clear, if 
not reversed, in research with PD. Specifically, higher education 
[4,9] and higher-status occupations [4—6] have been associated 
with a higher risk of PD or the results were null [10]. These 
findings are not well understood. Both education and occupational 
status may differentiate individuals based on intellectual engage- 
ment in various aspects of life. Occupational complexity is 
another, more refined way to measure intellectual engagement 
at work. 
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The concept of occupational complexity [1 1] captures cognitive 
effort required to perform a certain job, measured along the 
dimensions of complexity of work with data, people, and things. 
High occupational complexity is sometimes (but not always) 
reflected in higher education and/or higher occupational status. 
Building on research pointing to the association between high 
education level, high occupational status and greater risk of PD, 
we investigated the role of occupational complexity in risk of PD. 
Further, we used the unique opportunity to use data from the 
Swedish Twin Registry to assess directly whether any observed 
findings may be explained by familial factors shared by twins. 

We hypothesized that higher occupational complexity would be 
associated with increased risk of PD. We expected these results to 
exist for complexity of work with data and people in particular, 
which are often found to relate to health outcomes more strongly 
than complexity with things. Finally, we expected that any results 
would be attenuated in the twin pair analyses compared to the 
analyses with the entire sample, suggesting that familial factors 
such as inherited traits contributing to occupational selection or 
other factors shared by twins potentially underlie the results. 

Materials and Methods 

Study population 

The study is in compliance with the Swedish Act concerning the 
Ethical Review of Research Involving Humans. The study was 
approved by the regional Ethics Committee in Stockholm. Patient 
information was de-identified prior to analysis. As this was a 
register-based study, it was not possible to obtain written consent 
for all study participants. 

The Swedish Twin Registry (STR) [12] is the largest 
population-based twin registry in the world, containing informa- 
tion on Swedish twins born 1886 and onwards. In the 1960's, 
same-sexed twins born 1886-1925 were contacted with several 
questionnaires regarding demographic, medical and life-style 
factors. In 1973, same-sexed twins born 1926-1958 were 
contacted with a similar questionnaire. 

Between 1965 and 1990, national censuses were carried out 
every fifth year by Statistics Sweden, collecting information on 
current occupation, income, and type and size of habitation. 
These censuses have been linked to the STR through the unique 
10-digit personal identification number assigned to all Swedish 
residents. 

For the current study, we considered same-sexed twins of the 
STR born 1886-1950. Information from the 1970 and 1980 
censuses were used for those born between 1886-1925 and 1926- 
1950, respectively. Twins born 1886-1925 were included in the 
study population if they completed the first questionnaire in 1 96 1 
and were alive and without PD, according to the National Patient 
Register (NPR) and the Cause of Death Register (CDR), during 
the 1970 census. Twins born 1926-1950 were included if they 
completed the questionnaire in 1973, and were at least 30 years 
old, alive and without PD, according to the National Patient 
Register (NPR) and the Cause of Death Register (CDR), during 
the 1980 census. 

The original cohort contained 49,585 individuals. Of these, 
2,925 individuals were excluded because they died prior to the 
study period, were lost to follow-up, had PD before baseline, or it 
was not possible to link their data. Further, 17,882 individuals 
were excluded due to lack of complete data on the occupation or 
smoking variables. Of the people with missing occupational data, 
75% were women. Therefore, sex was included as a covariate in 
our analyses. The primary reason for lack of smoking data was 
failure of the participants to return a completed survey. The final 



study population consisted of 28,778 twins with 10,899 complete 
twin pairs. 

Occupational Complexity 

Information on occupation was obtained from the question on 
current occupation in the 1970 or 1980 Swedish census. For twins 
born 1886-1925, occupational data were obtained from the 1970 
Census, when the respondents were between 45 and 82 years old 
(M = 54.3 years, SD = 6.1). For twins born 1926-1950 occupa- 
tional data were obtained from the 1980 Census, when the 
respondents were between 30 and 55 years old (M = 41.3 years, 
SD = 7.1). Occupations were coded on a 3-digit level, based on the 
Nordic Standard Occupational Classification of 1965 ('Nordisk 
yrkesklassificering'; NYK), adapted from the International Stan- 
dard Classifications [13]. 

A job exposure matrix(JEM) was used to assess occupational 
complexity. Occupational complexity scores reflect those created 
in the fourth edition of the Dictionary of Occupational Titles [14] 
derived through extensive site observation of over 12,000 
occupations by qualified job analysts at multiple sites throughout 
the US [14]. Each occupation was rated regarding complexity of 
work with data, people, and things. See Table 1 for a description 
of the characteristics of work performed at each level of 
complexity. 

A study using a large sample of the US labor force in 1971 was 
conducted to assess inter-rater agreement of the initial ratings 
across occupations [15]. Kappa reliability estimates were 0.85 for 
complexity of work with data, 0.87 for complexity of work with 
people, and 0.46 for complexity of work with things [15]. Scores 
range from 0.0 to 6.0 for complexity of work with data, 0.0 to 8.0 
for complexity of work with people, and 0.0 to 7.0 for complexity 
of work with things, with lower scores reflecting higher complexity. 
To make the results more easily interpretable, we reverse scored 
the complexity variables so a higher score reflects higher 
complexity (Table 1). For example, physicists are relatively high 
in complexity with data, while lawyers are high in complexity with 
people, and precision tool makers are high in complexity with 
things. In order to address the positively skewed distribution of 
complexity of work with people, this variable was log transformed. 

To enable application of the JEM to Swedish occupational data, 
two independent raters, one in Sweden and the other in the 
United States, matched all the occupational categories from the 
1970 and 1980 Swedish Census to corresponding categories from 
the 1970 US Census. The raters initially agreed on 90% of the 
matches and the differences were discussed by the raters until they 
were resolved. Details regarding this procedure are published 
elsewhere [8]. 

Smoking Status 

Due to previous research showing the importance of smoking as 
a protective factor in PD [9], smoking status was included as a 
covariate to control for potential confounding. Information was 
taken from the questionnaires sent in 1967 and 1970 for twins 
born 1886-1925 and in 1973 for twins born 1926-1950. Smoking 
status was categorized as smoker vs. not smoker at the time of the 
questionnaire. 

Parkinson's disease ascertainment 

Incident cases of PD were identified through cross-linkage 
between the STR, the National Patient Register [16] (NPR) and 
the Cause of Death Register [17] (CDR). The NPR was initiated 
by the National Board of Health and Welfare in 1 964 and contains 
information on hospital discharge diagnoses from all hospitals in 
Sweden (mandatory reporting for all public and private hospitals). 
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Table 1. Description of the type of task associated with each level of complexity. 





Score range 


DATA 


Score range 


PEOPLE 


Score range 


THINGS 


5.6-6.0 


Synthesizing 


7.6-8.0 


Mentoring 


6.6-7.0 


Setting up 


4.6-5.5 


Coordinating 


6.6-7.5 


Negotiating 


5.6-6.5 


Precision Working 


3.6-4.5 


Analyzing 


5.6-6.5 


Instructing 


4.6-5.5 


Operating 


2.6-3.5 


Compiling 


4.6-5.5 


Supervising 


3.6-4.5 


Driving 


1.6-2.5 


Computing 


3.6-4.5 


Diverting 


2.6-3.5 


Manipulating 


0.6-1.5 


Copying 


2.6-3.5 


Persuading 


1 .6-2.5 


Tending 


0.0-0.5 


Comparing 


1 .6-2.5 


Speaking/Signaling 


0.6-1.5 


Feeding/Offbearing 






0.6-1.5 


Serving 


0.0-0.5 


Handling 






0.0-0.5 


Taking instructions 







Note. Shows the reverse scored range for the complexity levels and descriptions. Adapted from Ross and colleagues 2005 [8], 
doi:1 0.1 371 /journal.pone.01 06676.t001 



The register was gradually expanded and full national coverage for 
somatic disorders was obtained in 1987 [18]. Since 2001, the NPR 
also includes outpatient consultations with specialist physicians at 
hospital-based clinics. For the present study, data were available 
through 2009. The CDR contains information on all deaths of 
Swedish residents since 1 96 1 , including information on underlying 
and contributory causes of death. Causes of death data were 
available through 2008 and information on vital status through 
2009. Diagnoses were coded according to the International 
Classification of Diseases (ICD) in both the NPR and the CDR. 

Cases were defined as PD if there was a PD diagnosis in the 
NPR or the CDR. Both primary and secondary diagnoses, as well 
as underlying and contributory causes of death, were considered. 
Date of diagnosis was defined as the first date of diagnosis in the 
NPR, or the date of death in the CDR for those cases only 
identified at death. ICD codes used to define PD were: 342.00 
(ICD-8); 332.0 (ICD-9); G20 (ICD- 10). PD diagnoses in the 
inpatient register and the CDR have been validated against 
diagnoses based on clinical examinations: The positive predictive 
value(PPV) for PD diagnoses was 70.8% in the in-patient register 
and 66.7% in the CDR; the sensitivity was 83% in the two 
combined [19]. 

Statistical Methods 

The study population was followed from January 1 st 1971, for 
those born 1886-1925 or January 1 st 1981, for those born 1926- 
1950, until ascertainment of disease diagnosis, death, or December 
31 st 2009, whichever came first. Incidence rates (IR) of PD were 
calculated as events divided by time-at-risk, per 100,000 person 
years(PYR) with 95% confidence intervals (CI) assuming a Poisson 
distribution. 

We used two approaches to data analysis. In the first approach 
(Model 1) the association between occupational complexity and 
PD was modeled using Cox proportional hazards regression 
yielding hazard ratios (HR) with a 95% CI. Robust standard errors 
were used to adjust for any within-pair correlation and age was 
used as the underlying time scale in order to account for both age 
and follow-up interval within the model. We included sex and 
smoking as covariates in the model. We analyzed all individuals, as 
well as men and women separately. Occupational complexity was 
analyzed as a continuous variable. 

In the second approach (Model 2) we conducted a twin pair 
analysis using a stratified Cox proportional hazards regression, in 
which each stratum comprised a twin pair. This method allows for 
familial factors (genetic and early-environmental) shared by twins 



to be controlled for as part of the design. Thus, if the results are 
attenuated compared to the un-stratified design it would suggest 
familial factors account for part of the association. 

Finally, we conducted sensitivity analyses to examine differences 
by cohort, education, length of time in study and between people 
who were diagnosed while alive vs at death. Data analyses were 
performed with SAS software, Version 9.2 of the SAS System for 
Windows [20]. 

Results 

Demographic characteristics can be found in Table 2. We 
identified 433 incident cases of PD, 262 men and 171 women. The 
overall incidence rate of PD for the cohort was 56.0/100,000 
person/years (PYR; 95% CI 51.0—61.5) which is similar to other 
studies [21]. Incidence rates were nearly undetectable until age 65 
when, as expected, they started to rise exponentially and reached a 
level of 135.9/100,000 PYR in the oldest age group (90+ years), 
which is somewhat lower than other European studies [21]. 

In the first approach (see Table 3, Model 1), analyses indicated 
an association between higher occupational complexity of work 
with data and higher risk of PD, such that each additional level of 
complexity was related to a 7 % increase in risk, or an increase in 
risk by 42% between the lowest and the highest level of complexity 
with data. There was a trend suggesting higher occupational 
complexity with people was associated with a 10% increase in risk 
of PD but this was not significant (p = .05). When the analyses were 
stratified by sex, occupational complexity with data and people 
were risk factors in men but not women (p's>.05). No significant 
relationships were found between occupational complexity with 
things and risk of PD (p's>.05). 

In the analyses stratified by twin pair (see Table 3, Model 2), the 
associations between occupational complexity with data and 
people were no longer significant. However, the HR for 
occupational complexity with data did not change in these 
analyses. Conversely, the HR for occupational complexity with 
people dropped from 1.07 to 1.03. When these analyses were 
further stratified by sex the HR for occupational complexity with 
data dropped from 1.07 in men in the unstratified analyses to 1.00, 
but increased from 1.06 to 1.25 in women. The HR for 
occupational complexity with people decreased in men from 
1.13 to 1.02, but increased in women from 1.03 to 1.06. However, 
the confidence intervals in the stratified and unstratified models 
overlapped. 
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Finally, we conducted several sensitivity analyses. We conducted 
analyses separately for twins born 1886-1925 and 1926-1958 and 
found results similar to the overall analyses (for twins born 1886- 
1925,: HR= 1.07, CI = 1.00-1.14, for complexity with data-all; 
HR = 1.09, CI= 1.00-1.18, for complexity with data-men only; 
HR= 1.18, CI = 1.03-1.34, complexity with people-men only; for 
twins born 1926-1958: HR=1.09, CI = 0.98-1.21, jb>0.05, for 
complexity with data-all; complexity with data-men only, 
HR=1.05, CI = 0.93-1.20, p>0.05; HR= 1.06, CI = 0.87-1.28, 
p>0.05 for complexity with people-men only. 

We also examined the effect of education (basic education vs. 
higher than basic) as a covariate, exploring education as a proxy 
for those more likely to come into contact with the health care 
system. The HRs were similar as those from the original analyses, 
although the confidence intervals widened, reducing the statistical 
significance of these results (See Table 3). 

We excluded people who were in the study less than 15 years to 
reduce possible bias by reverse causation (n = 67); complexity with 
data-aU, HR= 1.09, 95% CI = 1.03-1.16; complexity with data- 
men only, HR= 1.13, 95% CI = 1.06-1.20, and complexity with 
people-men only, HR=1.16, 95% CI = 1.02-1.33. We also 
excluded people who were solely identified through the CDR 
(n = 34) to eliminate any possible differential bias between 
individual diagnosed while alive and those diagnosed at death; 
complexity with data-aU, HR= 1.08, 95% CI = 1.02-1.14; com- 
plexity with data-men only, HR= 1.08, 95% CI= 1.01-1.16, and 
complexity with people-men only, HR= 1.07, 95% CI =1.00- 
1.15. 

Discussion 

Several studies have shown a relationship between education or 
high status occupations and higher risk of PD [4,9]. We extend this 
research by examining the association between occupational 
complexity and risk of Parkinson's disease, presuming that 
individual dimensions of occupational complexity, which reflect 
education and occupational status to some extent, may shed light 
on the previous results by examining cognitive effort required at 
work. We found that higher occupational complexity with data 
and people were associated with increased risk of PD which 
appeared particularly high for higher complexity of work with 
people among men. Occupational complexity with things was 
unrelated to PD, probably due to its low reliability [15] or its 
inability to capture occupational characteristics relevant to risk of 
neurodegenerative diseases such as dementia [8] and PD. 

Analyses stratified by twin pairs yielded non-significant associ- 
ations, although the HR estimates were not obviously attenuated. 
Therefore familial factors do not appear to explain our results. 
However, twin-stratified analyses restricted to men suggested 
potentially greater familial influence in men. Thus, given the 
difference between the non-stratified and stratified results, familial 
factors such as personality and other factors affecting occupation 
selection, can be presumed to partially underlie the reported 
findings. Personality figures prominendy within the occupational 
complexity conceptual framework [1 1] as a predisposing factor to 
seek occupations of certain type and level of complexity, therefore 
we do not assume a causative relationship between occupational 
complexity and risk of PD. 

Many of the occupations shown to increase the risk for 
developing PD by previous research [4], such as business, 
healthcare, teaching, religious or legal occupations, would be 
coded as high in occupational complexity with data and people. 
Although not entirely consistent, epidemiological evidence indi- 
cates a pattern in which certain pre-morbid personality traits such 



as being industrious, punctual, inflexible, ambitious, orderly, and 
risk avoidant are common among individuals who develop PD 
[22]. Occupations that involve high occupational complexity with 
data and people are likely to be sought by someone who meets at 
least some of these characteristics. Alterations in the dopaminergic 
pathway [23,24], whereby certain personality types exhibit 
dopamine deficiency, offers another potential underlying mecha- 
nism for our observed association between occupational complex- 
ity and PD. However, not all research supports the hypothesized 
relationship between PD and personality [25] . The interaction of 
personality and occupation in relation to PD risk deserves further 
attention. 

Occupations requiring higher complexity with data and people 
may also involve greater levels of certain types of stress, such as 
interpersonal conflict associated with mentoring [26] in jobs such 
as managerial positions or professors [27]. Elevation of stress levels 
can cause increased levels of glutamate [28], which has been 
implicated in PD [29]. 

Limitations of the study include that occupation was only 
measured at one time point instead of over time. However, 
substantial job mobility could not be expected in this cohort [30]. 
Additionally, register-based PD diagnoses do not identify incident 
cases of PD by strict definition; however, they are a good 
approximation and there were only 34 cases identified solely by 
the CDR, and removing them from the analyses did not change 
the results. Registers have a PPV of 70% and it is highly unlikely 
that any misdiagnosis is differential by complexity score. In 
addition, there could be concerns about reverse causation bias, in 
which having preclinical signs of the disease may have led to 
different occupational choices; however, our sensitivity analyses 
which excluded people who were in the study less than 15 years 
suggested that any bias was negligible. 

Another limitation is the potential for medical surveillance bias, 
in which individuals with higher education or higher complexity 
occupations may be more likely to seek medical care. As there is a 
universal health care system in Sweden with equal access to health 
care across the population, there is a reduced pro-rich bias in 
obtaining healthcare [3 1] . Further, we controlled for education as 
a proxy for those more likely to come into contact with the health 
care system, but these results were similar to the main results 
(although statistical significance was reduced). These results 
indicate that occupational complexity contributes to PD risk 
beyond the effects of education. Finally, PD is chronically 
progressive and causes quite notable and disabling symptoms, 
minimizing opportunity for surveillance bias. 

Strengths of the study include the combined use of question- 
naire data with comprehensive record of PD diagnoses from the 
national health registers, with the questionnaire allowing for the 
assessment of smoking and education. Further, using the STR 
allowed us to control for the effect of genetic and early life factors 
using the twin pair analysis. 

In conclusion, by investigating individual dimensions of 
occupational complexity as a risk factor for PD, we provide a 
novel way of exploring previously reported associations between 
education, occupation, and risk of PD. People with an occupation 
high in complexity of work with data or people may be more likely 
to develop PD, with the increase in risk about even for these 
dimensions, albeit somewhat higher for complexity of work with 
people in men. In addition, we provide some evidence that familial 
factors, such as personality and other factors affecting occupational 
selection, may play a role in these results, suggesting that the 
underlying mechanisms may be inherent, at least in men. Further 
work to uncover the possible mechanisms of this association may 
be imperative for the success of any future intervention. 
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